Violation of the Equivalence Principle in Modified Theories of Gravity 
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We study modified theories of gravity of the f{R) type in Palatini formalism. For a generic f(R) 
lagrangian, we show that the metric can be solved as the product of a scalar function times a rank- 
two tensor (or auxiliary metric). The scalar function is sensitive to the local energy-momentum 
density. The auxiliary metric satisfies a set of equations very similar to Einstein's equations and, 
for weak sources, it can be approximated by the Minkowski metric. According to this, the metric 
coupled to the matter strongly departs from the Minkowskian one in the neighborhood of any 
microscopic physical system. As a consequence, new gravitationally-induced interactions arise and 
lead to observable effects at microscopic and macroscopic scales. In particular, test body trajectories 
experience self-accelerations which depend on the internal structure and composition of the body. 
These facts make very unlikely the viability of Palatini f(R) models designed to change the late-time 
cosmic evolution. 

PACS numbers: 98.80.Es , 04.50.+h, 04.25. Nx 



Modified theories of gravity have been object of in- 
tense study in recent years due to their ability to predict 
cosmic acceleration [jj without the need for sources of 
dark energy. In this sense, special attention have re- 
ceived those theories of the f(R) type, where / is a non- 
linear function of the scalar curvature R, with nonlinear 
terms relevant at low (cosmic) curvatures. For a given 
lagrangian /(-R), such theories can be formulated in two 
inequivalent ways, namely, in metric and in Palatini for- 
malism. In the former case, the connection is denned 
as the usual Levi-Civita connection. In the latter, the 
connection is regarded as independent of the metric and, 
therefore, must be determined by solving its correspond- 
ing field equations. According to this, the metric formal- 
ism leads to a system of fourth-order partial differential 
equations for the metric, whereas the Palatini formalism 
leads to second-order equations. Only for the linear la- 
grangian, f(R) = R — 2A, the two formalisms lead to the 
same equations of motion, which are those of General 
Relativity (GR). 

Though much has been written about the properties of 
theories in cosmological applications (see (J [3[ and 



for metric and Palatini formalism, respectively), we 
do not yet have a clear understanding of their properties 
in other regimes. Let us focus on Palatini theories. It is 
easy to verify (see below for details) that, for arbitrary 
/(i?), Palatini theories have the same type of vacuum 
equations as General Relativity with a cosmological con- 
stant. In particular, for a spherically symmetric star like 
the sun, the metric outside the star can be written as a 
Schwarzschild-de Sitter solution 



dr 2 

ds 2 Sds = g^dx^dx" = -A{r)dt 2 + — — + r 2 dtt 2 



A(r) 



(1) 



with A(r) = 1 - 2GM Q /r - Ar 2 /3, where A, G and M Q 
represent a cosmological constant, Newton's constant 
and the mass of the star respectively. The well-known 
model f(R) = R — ^ / R of Carroll et al. @, like any 
other f(R) model, admits such solutions and leads to 



a tiny A ~ /i 2 . One is then tempted to conclude that 
this model, and any other model with a small A, is 
compatible with solar system observations 0,[1,0| since, 
for sufficiently small A, (fT]) is virtually undistinguishable 
from the Schwarzschild solution of GR, which passes 
all observational tests. However, though the vacuum 
solutions are very well known, the equations in the 
presence of matter have not been studied in detail yet. 
The relation of the integration constant Mq with the 
matter sources, for instance, is still unknown. And this is 
an important aspect, since Palatini theories only depart 
from GR in the presence of sources, as is well known in 
cosmology. Some results in this direction were found in 
the study of the Newtonian and/or post-Newtonian lim- 
its in [8, 9]. All those works agree, with small differences, 
in that in addition to the usual Newtonian potential 
M/r the Newtonian limit has a density-dependent term. 
Such term was considered as irrelevant and negligible 
in some of those works. However, it was remarked in 
@ that it could be potentially dangerous. The reason 
being that it generates accelerations which depend on 
the gradient of the local matter density and which 
do not decay with distance. Furthermore, when an 
atomic/microscopic description of matter is considered, 
that density-dependent term seems to break the pertur- 
bative approach itself, which is a disturbing property. 
These facts suggest that further work is needed in or- 
der to fully understand the dynamics of Palatini theories. 



In this work we study and solve the equations of mo- 
tion of Palatini f(R) theories in the presence of matter 
in some cases of interest, namely, spherically symmetric 
distributions and also for very weak sources without any 
particular symmetry. From our analysis it follows that 
strong gravitational effects arise at microscopic scales 
which arc in conflict with the theoretical foundations 
of these theories and with the experimental evidence 
supporting the Einstein equivalence principle and our 
knowledge of elementary-particle physics. These facts 
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strongly suggest that f(R) models designed to modify 
the gravitational interaction at low curvatures should be 
ruled out. 



make our notation clearer, since i M „ and <7 M „ are con- 
formally related, it follows that TZ(T) = g^R^T) and 
R{b) — 9^Rpv(g) are related by 



Let us begin by defining the action of Palatini theories 
1 



2k 2 



trxy/=jjf(R) + S m \g^M (2) 



Here f(R) is a function of R = g» v R pv (T), with R^T) 
given by 



pA pp 
1 pp L v\ 



1 vp L p,X 



(3) 



where r* is the connection. The matter action S m de- 
pends on the matter fields ip m , the metric g^ and its 
derivatives, but not on r^„. Note that is regarded 
as a field of the gravity sector and, therefore, cannot ap- 
pear in the matter sector, since only the metric is allowed 
to couple to matter according to the postulates of metric 
theories of gravity Varying ((2]) with respect to the 
metric we obtain 



1 



f(R)R flv (T) - -f{R)g^ = k 2 T^ v 



(4) 



where f'(R) = df/dR. Note that the trace of © 

f'(R)R-2f(R) = K 2 T, (5) 

implies an algebraic relation between R and the trace T. 
The solution to this algebraic equation will be denoted 
by R = 1Z(T), The variation of © with respect to 
must vanish independently of ^ and gives 



(6) 



where /' = f'(lZ(T)) is also a function of the matter 
terms. Using an auxiliary tensor t pv = fg^u, © can be 
readily solved [l2| . The solution states the compatibility 
between and the metric t pu . In other words, can 
be written as the Levi-Civita connection of 



^ pv r, {dptpv + dytpn 



&ptpv ) 



(7) 



Inserting this solution for , written in terms of g^ u 
and f'(R(T)), in gj we obtain 

1 



Rpv{g) - -g pv R{g) 



1 



2(/') 2 



(8) 



where R^i,(g) and R(g) are computed in terms of the 
Levi-Civita connection of the metric g M „, i.e., they rep- 
resent the usual Ricci tensor and scalar curvature. To 



K(T) = R(g) 



2/'(T) 



f'(T) 



d x f'(T)d x f(T) - 77 ^_D/'(T) 

(9) 

where, recall, f'(T) = f'(lZ(T)) is a function of T, which 
means that the f'(T) terms in (|SJ) act as additional mat- 
ter sources. The matter terms 3 M /'(T) and □/'(T) make 
difficult the analisys of to obtain g^ v . In this sense, we 
find very useful to shift the problem and find solutions 
for the conformally related metric t^ t „ = (f>gu V , where 
4> = /'(T)//'(Q) is a dimensionless factor which becomes 
unity in vacuum, T = 0. The equations of motion for g^ v 
boil down to 



Gl{g) 



2k 2 " 



(10) 



where f a = /'(0) and we have used the shorthand no- 
tation V(T) = Kf'(T) - f{K). The above equations 
of motion for are considerably simpler than those 
for g^ v , which makes our task easier. In particular, 
for spherically symmetric systems we can use the ansatz 
ds 2 



ds 2 



g^dx^dx" 
1 



Airy^dt 



1 g ^ v dx^ dx v with 
1 



Mr) 



dr 2 + r 2 dfl 2 



which leads to 



2d<$> 

r dr 

1 rf(r[l - A]) 
r 2 dr 



f^ 2 

K 2 



1 r ± t 



A 



f^ 2 



t; 



V(T) \ 
2n 2 J 



Defining now A(r) = 1 - 2GM(r)/r in |T 
rewrite M(r) and $(r) as 



M(r) 
$(r) 



„2 



dx x 



n 



V{T)/{2n 2 



r 


■ rpr T't' 


/ dx x 
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<j) 2 A 



(11) 

(12) 

(13) 
we can 

(14) 
(15) 



Outside of the sources, (T^ v = 0, <j> = 1), the above equa- 
tions can be readily integrated leading to 



M(r) = M + 
$(r) = $ 



V{0)r 3 



(16) 
(17) 



and we recover the Schwarzschild-de Sitter solution JT]), 
with 2 A = V(0)/ Jq and 4>o =constant can be eliminated 
by a convenient redefinition of the time coordinate. 
Note that the constant M must be fixed by matching 
the interior and exterior solutions. Now, since the 
interior solutions depend on the details of T£ and the 
particular f(R) lagrangian considered through <fr(T) and 
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V(T), the value of Mq (and also $o) will depend on 
the internal structure and composition of the object 
(for a recent example see [HI). I n other words, a given 
amount of matter-energy, J drr 2 Tf, can lead to different 
external gravitational fields. This contrasts with Gen- 
eral Relativity, where the external field generated by a 
spherically symmetric system only depends on the total 
matter-energy and not on the structure or composition 
(in GR 4> = 1) V(T) = 0). This behavior was already 
reported in Q within the post-Newtonian regime. 

We have just seen that outside of the sources (fTTj) 
turns into ([lj. Let us consider what happens if we 
place a test body with m <C Mq within the external 
field of Mq. This test body could be anything from a 
laboratory-sized object to a single atom. The metric 
nearby this body is given by |T]). However, the metric 
inside of the test body, where T ^ 0, is given by JTTj) 
with 4>(T) 7^ 1. In the 1/R model, for instance, we 
find 6 = 1— ; 1 where a = — TlT r and 

2[l + Vl + 12/<r 2 ]' ' 

T c = li 2 /k 2 = p^c with ~ 1CT 26 g/cm 3 . When T 
drops below T c (outside the body) we have <f> — ► 1. For 
T above T c (inside the body) we get <j) ^ 3/4. We 
then see that the line element (fTl"]) changes suddenly 
from ds 2 = ds 2 sds to ds 2 — |c?s| dS when going from 
the outside to the inside of the test body. Actually, this 
happens in general when going from the outside to the 
inside of atoms, which leads to strong oscillations in the 
metric at microscopic scales both within the test body 
and within the central object characterized by Mq. 
This unusual and highly non-perturbative behavior, 
which has no precedent in the literature on alternative 
theories of gravity, is due to the dependence of the 
metric on the factor <j>(T), which is a function of the 
local energy-momentum density T. Any theory in which 
the lagrangian f(R) be sensitive to some low curvature 
scale will lead to a <j>{T) characterized by two limiting 
values, (j> — > 1 for small T and 4> — > <p c for large T, as 
compared to the characteristic scale T c . Therefore, the 
strong jumps in the metric discussed here for the l/R 
model occur in all f(R) models aimed at changing the 
late-time cosmic expansion. 

To see in more detail the effects of the factor (f>(T), 
we will consider a portion of space-time containing the 
above test body. Assume that we can take coordinates 
in which the line element ([I]) becomes Minkowskian away 
from the body, i.e., ds 2 m T)ni/d£ ft d£ 1 ' ■ This choice of co- 
ordinates eliminates the external field generated by Mq. 
We could also have assumed no external field for sim- 
plicity. The metric g pv in the region close to our test 
body can then be computed using again the decomposi- 
tion g^ v — 4>g pi ,. For a sufficiently light test body, such 
as an atom, the metric g pv can be well approximated by 
?7 M „ . This is so because g^ is governed by (fT0|) , which is 
dynamically very similar to General Relativity. In fact, 
the term k 2 /(/q0 2 ) can be seen as an effective (density 
dependent) Newton's constant, and V(T) as some (den- 



sity dependent) cosmological constant. For lagrangians 
of the form f(R) = R + eg(R), with e some small param- 
eter, we find that V(T) ~ e. Thus V(T) is very small (of 
cosmological order) and can be neglected in the region 
of interest. Like in General Relativity, the contribution 
of T" is also negligible for weak sources. Therefore, (jTUJ) 
becomes G"(g) ~ and leads to g^ v ~ rj^. The metric 
generated by our test body (atom) is then 



flV(z) = 1 [T{x)]rj ftl/ 



(18) 



Note that the factor _1 [T(a;)] is due to the derivatives 
of f'(T) on the right hand side of ((5J) which, unlike 
and V(T), cannot be neglected. We thus see that the 
metric in the neighborhood of any small physical system, 
described by a certain probability density distribution, is 
not in general the Minkowski metric. This result is very 
important and leads to a number of inconsistencies and 
pathological effects, which must be added to those dis- 
cussed above. 

Let us first focus on the theoretical inconsistencies. The 
action @ was constructed according to the postulates 
of metric theories of gravity, namely, i) spacetime is en- 
dowed with a symmetric metric g pv , ii) the world lines of 
test bodies are geodesies of that metric, and hi) in local 
freely falling frames the nongravitational laws of physics 
are those of special relativity. These postulates tell us 
how to introduce matter in a theory of gravity based on 
geometry: we must take the Minkowskian theory and 
change r)^ by g^ and by V M , with V Q 5 M „ = 0. This 
prescription should guarantee that the non-gravitational 
laws of physics of Minkowski space were recovered in lo- 
cal freely falling frames. In our case, however, even in 
local frames in which external gravitational fields have 
been screened, we do not recover the Minkowskian met- 
ric, which violates postulate hi). 

On the other hand, since these postulates assume that 
the Einstein Equivalence Principle (EEP) is valid, it fol- 
lows that Palatini f(R) theories must be in conflict with 
the experimental evidence supporting the EEP (see [ll| 
for details). Let us thus focus on the weak equivalence 
principle (also known as universality of free fall, and a 
key piece of the EEP), which states that the trajectory 
of a freely falling test body in a local inertial frame is 
unaccelerated and independent of the internal structure 
and composition of the body. From (fTS)) we see that the 
geodesic equation 



d 2 e 

dT 2 



CUg) 



dj a d^ 
dr dr 



0. 



(19) 



where C^(g) = ^- {d a g p p + d f jg pa - d p g a p), leads to 



dv l 
~~dt 



-(1 - \v\ 2 ) [?M In + ^1110] (20) 



where v l = d^ l /dt. This equation shows that if d p (f>(T) ^ 
0, the test body will feel self-accelerations which, in addi- 
tion, will depend on its own internal structure and com- 
position. Now, the deep dependence of <p(T) on the mi- 
croscopic structure does not allow us to use an averaged 
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description of matter (perfect fluid) to study the behavior 
of (|20[) . One should use a microscopic description for the 
matter sources, as pointed out in [13| . Let us then con- 
sider the action of a free Dirac field in Minkowski space 
and use the EEP to construct the corresponding curved 
space theory to be put in ([2|). The result is 



-g [WV^ + imiH>] » (21) 



where = (<9 M — T^)ip, is the spin connection, 

and {7^,7"} = 2g tiv are the curved-space counterparts 
of Dirac's constant matrices {7 a ,7 6 } = 2r] ab . Though 
the Minkowskian theory should be recovered in a freely 
falling frame, where external gravitational fields have 
been screened, from (fl"8| and (|2Tj) we find 



(22) 



where T depends on all the sources present, including ip 
itself, and c^ 1 / 2 is due to [lq . It is quite evident 



that (|22|) is not the Minkowski space theory we started 
with and that new gravitationally-induced (first-order) 



interactions arise 17[ via the various 4>{T) terms in 



This modification of the microscopic interactions 
is what eventually leads to the macroscopic violations of 
the universality of free fall observed in (j2"0)) . Note that 
this effect is different in nature from the Nordtvedt effect 
[i~4j . which is due to the coupling of the gravitational 
self-energy of a massive body to external gravitational 
fields (second order gravity-gravity interaction) . On the 



other hand, the fact that 0(T), in models relevant for the 
late-time cosmology, is sensitive to low energy-density 
scales, makes it difficult (or perhaps impossible) an 
expansion of (|22[) in powers of T to quantitatively study 
the predictions of particular models. This means that 
the non-gravitational laws of physics of Minkowski space 
are strongly modified and might even lose their predic- 
tive power. This problem also affects the gravitational 
laws, since the metric is sensitive to the microscopic 
structure of T M „, and we must solve the microscopic 
matter field equations to get T^ v . This fact, by the 
way, invalidates all cosmological applications of Palatini 
theories considered so far 0, Q, in which an averaged 
description of matter (perfect fluid) has been assumed. 
In summary, we have shown that the gravitational sector 
in Palatini f(R) theories has dramatic effects on the mat- 
ter sector. When the gravitational interaction is turned 
on, the Minkowskian laws of physics are completely lost 
due to the sensitivity of the spacetime metric to the local 
energy-momentum densities. This fact, besides being a 
theoretical inconsistency, implies observable effects such 
as violations of the EEP and new interactions among 
elementary particles. Furthermore, if the theory is sen- 
sitive to very low density scales, it might be impossible 
to use the theory to extract quantitative predictions (see 
[TH for an example). All this suggests that Palatini f(R) 
theories, at least in their current form, should be ruled 
out. Only those models in which <p{T) is only sensitive 
to very high energy densities, above the currently acces- 
sible experimental limits, might have a chance to survive. 
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This result is independent of the choice of conformal 
frame for the gravitational action 



